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Abstract 

Through a specific example of two-body color-favored charm decay, Df 4''^'^} we 
illustrate how an effective and complex (unitarized) oi, denoted by a!('^^^ , may be defined 
such that it includes nonf actor ized, annihilation and inelastic final state interaction (fsi) 
effects. The procedure can be generalized to color-suppressed processes to define an 
effective, and complex a!^'^^^ . We determine |a^''^''^'^| and, where relevant, \(i2'^^^\ for 
D Ktt, Kp, K*TT, — > rjir'^, rj-n^, rip^ ,rj ^ and for D^tt^ and from 

the hadronic and semileptonic decay data. 



^Present address: Nuclear Physics Division, Bhabha Atomic Research Centre, Bombay 400085, India. 



1 Introduction 



We begin with some definitions relevant to the hadronic decays of charmed mesons. The 
effective Hamiltonian for Cabibbo-favored charmed decays is given by, 

K"{AC = AS = -1) = Gf {Ci{sc){ud) + C2{sd){uc)} , (1) 

where Gp = ^^csVud the brackets (sc) etc. represent color-singlet (V-A) hadronic cur- 
rents with appropriate flavors and Vcs etc are the Cabibbo-Kobayashi-Maskawa (CKM) mixing 
parameters. Gp is the Fermi Weak coupling constant. Gi and G2 are the Wilson coefficients 
^ for which we take the values, 

Ci = 1.26 ±0.04, C2 = -0.51 ±0.05 . (2) 

The central values are taken from [0]; the error assignments are ours. The parameters ai and 
02 are defined as follows, 

I3'l,2 = Cl,2 ± , (3) 

where is the number of colors. 

The formula corresponding to (1) and (2) for Cabibbo-favored bottom decays are, 

H^f^iAB = AG = -1) = ^V,,V:, [G,{cb) {{du) + (sc)} ± G2 {{cuKdb) + {cc){sb)}] , (4) 

and i, 

Ci = 1.12 ±0.01, C2 = -0.27 ±0.03 . (5) 

Again, the error assignments are ours. 

In phenomenology as practiced until recently, it was found that the choice ai^2 = (^1,2 
worked as a reasonable approximation in the factorization scheme for charmed decays. The 
most successful example of this was the decay D — > Kn 0, ^. This led to the belief that 
Nc ^ 00 was a good approximation in charmed decays. This idea, when carried over to 
hadronic B decays failed as theory wanted 02 to be negative while experiments left no 
doubt that in B decays 02 was positive. 

Recently, it was shown that in the factorization hypothesis all commonly used models of 
hadronic form factors had difficulty in explaining the polarization data P, |^, || in 5 — > tpK* de- 
cay. It was subsequently shown in Q that even a modest amount (~ 10%) of nonfactorized con- 
tribution made all form factor models consistent with the polarization data. The consequences 
of nonfactorization in charmed meson decays have recently been explored in |l^, |T^, 0, |l3l and 



in B decays in [in 



If we use A^c = 3, we get from (2) and (3) at the charm scale, 

ai = 1.09 ±0.04 , a2 = -0.09 ±0.05 , (6) 

and at the bottom scale, 

ai = 1.03 ±0.01 , 02 = 0.10 ±0.03 , (7) 
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In 12 and 13 



it was shown that the inclusion of nonfactorized contributions allows us to 
define effective ai and 02 and that even a modest nonfactorized contribution in color- suppressed 
charm decays could lead to a^^ ~ —0.5, a circumstance mitigated by the large value of the 
ratio It has not yet been explicitly shown as to what effects are included in a^^ and : 
Are they complex, and if so, what makes them so? Where do annihilation processes fit in? 
What role do final state interaction (fsi) play? 



2 An Illustrative Example. 

We answer all the questions posed at the end of the preceding section with a specific example 
from charm decays: —>■ ^vr"*" . The reason for choosing this Cabibbo-favored decay is that 
it has only one isospin which makes the fsi calculation somewhat simpler. Before embarking on 
the details, let us introduce the following color Fierz identities, 

1 1 ^ 

isc)iud) = —{sd)iuc) + 2 E (sA"d)(MA'^c) 

c a=l 
1 1 ^ 

where A"^ are the Gell-Mann matrices. We adopt the following short-hand notations for the 
color-octet current products on the right hand side of (8), 

^if ^ = \Y. (^"^"c) (^>^''d) , ^(f ) = \y. ('^^""d) (uX'^c) . (9) 

Using (1) and (8), the decay amplitude (before fsi effects are brought into play) for 07r+ 
is given by, 

A{Dt - 07r+) = Gf {a, {<l>n+\{sc){ud)\Di) + C2 (07r+|ij(f) . (10) 

While the matrix element of H^^'^ is completely nonfactorized, the first term in (10) includes 
a (i) factorized (spectator) term, (ii) any nonfactorized contributions in addition to the fac- 
torized (spectator) term and (iii) a W-annihilation term which in turn has a factorized and a 
nonfactorized part. These individual contributions to the decay amplitude are parametrized as 
follows: 

[<Prr^\{sc){ud)\Dty'"'' = M2m^)e.pnX"^{0 , 
[<PiT+\{-sc){ud)\Dty^ ^ U2m^)e.pDA^''^ , 

, , \ ann 

[4>n^\{sc){ud)\Dt) ^ fnA2m^)e.pDAr , 

^^^\HL'^\Dt) = U2m,)e.pnA^''^ • (11) 

Here we have used the form factor notation of while A'q^"^ and Aq^^"'^ were introduced in 



r^. The superscript 'ann' stands for annihilation and e is the polarization four- vector for the 



'). Putting it all together, one can define an effective oi as follows, 

AiDj ^ 07r+) = GFalf^U2m^)E.pnX"^irnl) , (12) 
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where 



ann 



a\^^ =aAl + -rn; + + } ■ (13) 

Up to this stage, all quantities are taken to be real, including Aq™. Complex amplitudes will 
emerge as the result of fsi at the hadronic level. 

Consider now the final state interactions. For illustrative purposes we consider a two- 
channel model which is adequate to illustrate our ideas. Consider an inelastic coupling of 07r 
channel with G-parity even, to the G-parity even eigenstate of K^K*~^ and K*^K~^. Channel 
(j)7i will couple, among others, to rjp and rj'p channels also. Our intention is not to calculate 
numerically the effect of these channels but to illustrate how fsi enter the description. Both of 
these decays, — > K^K*^ and — > K*'^K^, are color-suppressed. Following an analogous 
procedure to the one that led us to (12), we find, 

A{Di ^ K'K*+) = GFal^^fK{2mK^)e.pnAo''''iO , (14) 

where 



ann 



TDa 

Aq [m^) a2 Aq (m^j 02 Jk Aq (m^j 

Here B^^^""^ , 5^^^"-^ and 5^"" are the analogues of aJ,^^"-^ and Ag'"" of (11). Similarly, 

A{Di ^ K*'K^) = Gpalff fK^{2mK^)e.pnF^^''{ml,) (16) 

where 



lann 



TDar 



The hatted quantities here refer to the decay channel K*^K^ . 
Now the eigenstates of G-parity are |[T5[] , 



\K*K)s,A = i {\K*+K^) ± \K+K*^)} 
v2 ^ ^ 



where the symmetric (antisymmetric) combination has G-parity even (odd). Thus it is only 
\K*K)s that couples to (pir. We note here one further point regarding the annihilation term 
in (15) and (17). The factorized part of the annihilation term in (15) is (^K^K*~^\{ud)\o'^ 
{0\{sc)\Df) which is proportional to the matrix element of the divergence of the axial part 
of (ud) current. Now, if the hadronic weak currents are only of the first class kind then the 
axial current has G-parity odd. As the sjTiimetric state, \K*K)s, has G-parity even, it requires 
that the factorized part of IJq"" in (17) be equal in magnitude and opposite in sign to the 
factorized part of Bq"-^ in (15). However, nonfactorized annihilation terms (for example, when 
the intermediate state in the direct channel is not a hadronic vacuum but a multigluonic state 



iTql) will frustrate this argument. 
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We, next, set up a coupled channel fsi between the decays Z)+ — >• 07r+ and {K*K)s 
following the formalism described in [|T^. Though the method of unitarization, the K- matrix 
method which amounts to retaining only the on-shell contribution from re-scattering loops, is 
not unique, it serves adequately to describe our ideas. 

We simplify our notations further by using the following short-hand notations for the thus 
far real amplitudes, 

A{Dt ^ 07r+) = s.pnA'^'' , (19) 

with 

A^^ = G^af ^^(2m^)<»^(m2) , (20) 

and 

A{Dt - iK*K)s) ^ e.pnA''''' , (21) 

where 

A^'^ = G^^^ {a," fKA^^^\ml) + aj^ Ik^F,^^^ (mj,,)} . (22) 

The two amplitudes, (19) and (21), couple via fsi and get unitarized. The unitarized decay 
amplitudes, , are given by [ |T7[| , 

= {l-ik^K)'^^ A , (23) 

where A is a column with entries A"^^ and A^*^, k'^ is a diagonal matrix with entries kf and 
/c|, ki and k2 being the center of mass momenta in the channels (pn and K*K respectively and 
K is the symmetric, real (2x2) K-matrix, 

^-Ce)- 

where a, b and c are assumed to be constants with dimensions GeV^'^. Note the appearance of 
fc^ as the appropriate threshold factor for P- waves in (23). 

The parameters of the K-matrix could be evaluated from the knowledge of the coupled 
channel scattering problem. In absence of this information, they remain undetermined in our 
case. Though, for our purposes the knowledge of the numerical values of the K-matrix is not 
necessary, we have ventured an estimate of the elements of the K-matrix later. 

Carrying out the unitarization of the decay amplitude as indicated in (23), we obtain a 
unitarized A^''^''^ which is complex and given by, 

A^iDj ^ 07r+) = G^af'^^^Ae.pB,(2m<^)<=^(m^) , (25) 



where 

= ^ U - ^klc + r^kih \ + ^^^Wr^ ) \ > (26) 
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with A = det(l - ik^K). 

If the fsi were elastic, b = c = and A = 1 — ikfa, we would have obtained 

a^" = -j^—e'' , (27) 



where 6 = tan~-'^(aA;J) is the elastic P-wave 07r scattering phase. 

Though numerical calculations is not the intent of this paper, a rough estimate of some 
of the K-matrix elements can be obtained in the manner done in [T^. For example, the off- 
diagonal T-matrix element,Ti2, representing the inelastic process 07r {K^K*^ + K*^K~^), in 
the K-exchange approximation, is given by, 

Tu{W, e) = 2^^pp ^^■P'^^K' -P. ^28) 
[Pk* -PnY - mj^ 

where W = m£,^ is the center of mass energy and we have used an SU(3)-symmetric Vector- 
Pseudoscalar-Pseudoscalar (VPP) coupling gypp given by [1^ 

4^ ^ 3.3 . (29) 

4:71 

The fact that the decaying particle has spin zero imposes simplifying helicity constraints 
on the vector particles in the process (pn —>■ {K*K)s- As the orbital angular momentum, L, is 
orthogonal to the scattering plane, it cannot have a component in the plane of scattering. This 
forces the helicities of and K* to be zero in the re-scattering process. We can then project 
out J = amplitude with = 0, Xk* = from (28) by using 

{Tu{W, ^))a„a,. = E (2^ + 1)<,A,. {0)%,x,, {W) (30) 
J 

Projecting (T'{^"°)oo from (28), we obtain 

(^i2~°)oo = ^^^-P {{EkEt, + -Ek*Es)Qq{z) 
m^niK* 3 

- — {p'^E^E^ + p'EkEk.)Qi{z) (31) 
pp' 

2 

+ -E^Ek.Q2{z)} 

where p and p' are the magnitudes of 3- momenta in the center of mass of (/)7r and K*K systems 
aX W = m£,^ and Qi{z)aie the Legendre functions of the Second Kind with the argument z 
given by 

z = ^^{2EkA + ml^ - ml) (32) 

Finally , we relate (T{^^'')oo to the off-diagonal K-matrix, through 

(T^it°)oo = 87riyA;iKi2fc2 (33) 



5 



where ki and ^2 are the eigen-momenta in the two channels and appear due to the P-wave 
nature of scattering in J = state. Numerically we obtain 

iT,i=')oo = l.58g'ypp, Ku = b = 2.73GeV-' (34) 

By a similar technique, one can calculate K22 through tt exchange in the diagonal channel 
{K*K)s {K*K)s which goes via K*+K° K*°K+ . We obtain 

(T/2=°)oo = 0.99(7^pp, K22 = c= l.TSGe^-^ (35) 

These parameters, b and c, are quite sizable. Their effect in the unitarization appears 
in dimensionless products of form k\h^ Zcffc and A;|c which are numerically 0.98, 0.87 and 0.57 
respectively. It is harder to estimate Kii{= a) as the elastic 07r scattering is an OZI- violating 
process. Thus, though our model is crude, it appears very likely that inelastic fsi could play an 
important role in (pn decay. 

A similar expression to (26) can be written down for Df ^ K*+K° (and for D+ K+K*^) 
which would define a^'^^^ and a^'^'^'^. 

One should view (25) as the defining equation for which includes all conceivable 

physics, is pro cess- dependent and complex. If we view a(''^^^ and a!^''^^^ in the manner we are 
proposing, then the comparison of two body hadronic decays of D and B mesons with semilep- 
tonic decays which in past has been claimed [0, |, |2y] to be tests of factorization, becomes 
merely determinations of |af''^'^'^| and \ci2'^^^\- 



3 Estimates of \a('^^^\ and |a^'^^^| from charm and beauty 
decay data . 

The procedure we have outlined above can be used in defining a(''^^^ and a^'^^^ in, say, Kti 
decays. There is an added complication here, that of two isospins in the final state. The fsi 
unitarization has to be carried out in each of the two isospin states separately. Nevertheless, 
one can define, following the same procedure as we have used for the simpler case of — » ^tt"*", 

A{D' ^ K%') = ^\a'^'^ff\f^{ml - ml)F^- {ml)e^^^^ (36) 

A{D+ K\+) = AiD"^ K-n+) + V2A{D^ i?°7r°) . 

In principle, |ai^'^"'^''^| can be determined by relating T{D'^ K~tc^) to r(Z}° K~l^v) and 
\cL2''^^^\ by relating T{D^ A'^vr") to T{D^ — > t:~1^v). Finally, 0oo is, in principle, obtain- 

able from T{D^ ^°7r"'").We determined the products \a('^^^\F^K{m^) and \ci2'^^ ^ \F^'^ {ra\) 
and the relative phase (0+_— 0oo) from the branching ratios B{D^ — > K^-k^), B ( 7r° ) 

and B{D+ ir°7r+) ^ with the result: 



D^Kn: \a^'^f^\F^i<{ml) = 0.767±0.014 

\(h'^"\F,^^irnj^) = 0.593 ± 0.038 (37) 
cos(0+_ - 0oo) = -0.867 ± 0.089. 
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If we use the experimental determinations |21| of -Fq^^(O) and -Fq^'^(O) from semileptonic 
decays (assuming monopole extrapolation), 



■^DK 



(0) 

Fo^-(0)/Fo^^(0) 



0.75 ±0.02 ±0.02 
1.0i[!;^ ± 0.04 
1.3 ±0.2 ±0.1 



[21] 

m 



(3^ 



23 



we obtain, 



1^2 I 



1.02 ±0.04 
(0.76t°:?^), 



(0.58 ±0.08) 



(39) 



In a^'^^^ above, the two values correspond to the two values of the form factor ratio 
Fq'^{Q) / Fq^{Q) given in (38) respectively. We have used a monopole extrapolation with pole 
mass 2.47 GeV in calculating F^'^{m\). 

A similar analysis of the branching ratios in D 



\a 



U,eff\ 

1 I 

U,eff 



l«2 
COS 



Fl 



K) 



Kp and i^*'7rleads to: 

1.097 ±0.069 
0.672 ±0.055 
-1.046 ±0.205 



(40) 



and 



D K*7r 



\a 







\0'2 
COS 



1 



{ml)= 1.138 ±0.070 
; ) = 0.747 ±0.061 
()oo) = -0.926 ±0.166 



m 



(41) 



From the form factors at = listed in we can calculate all the form factors needed by 
using monopole extrapolation for all of them except AQ^{m\) for which we adopt the theoretical 
value given in P| . The resulting a!(''^^^ and a!^''^^^ are : 



Kp: 



1.27 ±0.09 
0.93 ±0.08 



(42) 



and 



D ^ K* 



TX 



\0'2 I 



1.76 ±0.23 

(0.8tEJ:?^), (0.61 ±0.09) 



(43) 



The two values of a!^'^^^ correspond to the two values of the ratio Fq'"{Q)/ Fq^{^) respectively, 
given in (38). 

We end with a determination of the pro cess- dependent in — > t^tt^, r/V^, T^p"*" 

and r]' p'^ and — > D'n^ and D~p^ from hadronic and semileptonic data. For decays, 
we provide a calculation for rjir^ and rip~^ to illustrate the method, details of which may 

be found in Refs. and p4 . 
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The defining equation for a!('^^^ in —>■ rj7T~^ and rjp'^ decay amplitudes is obtained by 
simply replacing ai in the expression for the factorized amplitude by a(''^^^ as in (25). Thus 

A{Dt - w^) = GpC,{a1^^ff),^.U{ml^ - ml)F^'\ml) , (44) 
where F^"^ is the relevant form factor [Q, and in terms of flavor singlet-octet mixing angle 9p, 

Cr, = ( cos^p + ^ sin^p ) . (45) 



V2 



The resulting decay rate is, 



TiDj^Vn^) = {CM<.-rnl))\{ml^,mlml)\F['^^0)\' , (46) 



where A(x, y, z) = (x^ + y'^ + z'^ — 2xy — 2xz — 2yzY^'^ and we have used F^'''^{ml) ^ Fo^"''(0) = 
Fi^=''(0). Similarly, 

A{Dt ^ r^p+) = GFCr,{a';''"),,^{2mMe*.poF['''{ml) , (47) 



which leads to 



where n = 1(2) for a monopole (dipole) extrapolation of the form factor. Ai is taken to be 2.11 
GeV, the mass. 

Consider now the semileptonic decay rate for —>■ rjc^v pO|, which can be written as, 
T{Dt - ve^u) = ^^^Cl\F^''m'^Tln{mn^.m,,K,) , (49) 



where 



In{mn.,m„A,) = dq' {q^-Air ' ^^^^ 

j'From (46), (48) and (49), we can construct the ratios the ratios ^^^l and + ^^p^) 

)^^+| and \iai'''^^)^p 



involving the unknowns |(af''^'^''^)n7r+ 1 and |(af''^'^'^)np+ 1- By equating these theoretical ratios to 



the experimental ones, we can determine |af'^"'^''^| in these two decays. A similar method can be 
applied to decays involving rj' in the final state. 

We now turn to the experimental results we have used. Recently, CLEO collaboration has 
measured ||2^ the following ratios, 

T(Dt ne+u) 

-^-^ ^—4 = 1.24 ± 0.12 ± 0.15 , 

Y(D+ ^ 0e+z/) 

(51) 

T(Dt ie+u) 

— '—-4 = 0.43 ± 0.11 ± 0.07 . 

r{D+ ^ <Pe+u) 



If we combine this with the following measured ratios, 

r{Dt ^ 07r+) 



and 



r(^+ - 




nDt - 






nDt - 


r]'n+) 
















- Vp+) . 


nDt 





0.54 d= 0.05 d= 0.04 [21], (52) 



0.54 ± 0.09 ± 0.06 [26] , (53) 



1.2 ±0.15 ±0.11 [26], (54) 



2.86±0.38l^J^ [27], (55) 



3.44 ± 0.62i°:^^ [27] . (56) 



We obtain the following experimental ratios: 



and 



nDt - 




nDt - 


77V+) 


nDt - 


^ 77'e+i/) 


nDt - 


77P+) 


nDt - 




nDt - 





r{Dt r]'e+v) 



= 0.81 ± 0.23 , (57) 



= 5.17 ± 1.86, (58) 



4.27 ±1.13, (59) 



= 14.81 ± 5.81 , (60) 



The errors here are probably overestimated as we propagated all errors as if they were inde- 
pendent while some systematic errors in the products of ratios would cancel. 

Confronting the theoretical ratios to the experimental ones shown in (57) - (60) we have 
evaluated the following (we have used Kd = 0.975, U = 130.7 MeV and fp = 216.0 MeV ): 



|(af''^0»7vr+| = 0.89 ±0.13 (n = 1) 

= 1.08 ±0.15 (n = 2) , (61) 

(af''^0»7'7r+| = 1.56 ±0.28 (n = 1) 

= 1.68 ± 0.30 (n = 2) , (62) 
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\{ai'^")vP+\ = 1-49 ±0.20 {n = 1) 

= 1.55 ±0.20 {n = 2) , (63) 

\i(h'^^^)v'P+\ = 2.77 ±0.55 (n = 1) 

= 2.60 ±0.51 (n = 2) . (64) 

In most cases |af''^'^''^| has risen for a dipole form factor compared to the monopole, except 
for —>■ i]'p^. The reason is that the hadronic rate, T{D^ —>■ rj'p'^), rises more than the 
semileptonic rate, T{Df — > rj'e^v), when one goes from the monopole to a dipole form factor 
in contrast to the other cases. 

For D^TT^ and D~ p~^, we define ai'^^^ via, 

A{B' D-n^) = ^V;;K. K '^'Oz.. K - ' (65) 

A{B' D-p^) = '^^Oi,^ (2m,/,) e.^B^rK) ' (66) 

where we have used F^^{m1) ~ F^^{0) = Ff^{0). The hadronic rates are calculated from 
the above two equations and the semileptonic rate from an analogous formula to (49). For the 
experimental branching ratios, we used Ref. and evaluated |af''^'^'^| for B Dtx and Dp 
decays. We used four different extrapolations for the form factor Ff^^{q^): (i) monopole and 
(ii) dipole with pole mass 6.34 GeV 0, (iii) an exponential form as in Ref. [pH] , 

^(t = q^) oc exp {0.025(t - t^)} (67) 



where tm = {niB — in GeV ^, and (iv) a form advocated in 

2 



Ff^(t = g2) 0, exp \-(3- — ^1 (68) 



where (3 = 2p^ - I with p = 1.19 || and, 

_ / (m| + ml- q^) 
y = v.v = (69) 

The resulting |af'^'^'^| for each of the four form factor extrapolations are listed below. 
For B Dn: 

l^f^.e/Zj _ Q 92 ± 0.05 (monopole) 

= 1.0 ±0.06 (dipole) (70) 
= 0.91 ± 0.05 (using (67)) 
= 1.0 ± 0.06 (using (68)) 



For B Dp: 



\a^'^^^\ = 0.91 ±0.08 (monopole) 

= 0.99 ±0.09 (dipole) (71) 

= 0.91 ± 0.08 (using (67)) 

= 0.99 ± 0.09 (using (68)) 
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4 Summary 



In summary, the effective, and unitarized ai and a2 which are defined by the following pre- 
scription: The true decay amplitude is given by replacing oi and 02 by a'^''^^-^ and a^'^-'^'^ re- 
spectively in the factorized (spectator) amplitude. Defined in this manner, as we have shown 
systematically how these effective parameters get contributions from annihilation and nonfac- 
torizable processes as well as the final state interactions. As these effective parameters are 
process-dependent, the purported test of factorization that compares the hadronic rate to the 
semileptonic should be used, instead, as a tool to determine the modulus of these effective 
parameters. 

We determined |af''^'^'^| and |a^''^'^'^| in D — > Kn, Kpand K*t: decays using experimental 
input on formfactors ( with monopole extrapolation) as much as possible. The values of these 
parameters, particularly \a!^''^^^\, imply large departures from the factorization expectation 
when compared with ai and 02 given by (6) with Nc = 'i. 

For "^T^^ 1 and t]' ^ a{'^^^ was found to be significantly different from a\ of 

eq.(6) signifying that the simple factorization prescription would not apply to these cases. For 

— »• D~n~^ and D~p'^, however, we found |af''^'^''^| to be not much different from ai of eq.(7), 
especially for the form factor extrapolations given by a dipole and eq.(68), signifying that effects 
such as annihilation, nonfactorization and fsi play a less significant role in hadronic B decays. 

Finally, as emphasized in Refs. and 0], effective ai and 02 can be defined only for those 



decays whose amplitudes involve a single Lorentz scalar structure. Thus they can not be defined 
for decays of D and B mesons involving two vector particles in the final state. Consequently, 
our analysis applies only to those cases where the decay amplitudes involve a single Lorentz 
scalar structure. 
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